Mutually unbiased bases in dimension six: The four most distant bases 
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We consider the average distance between four bases in dimension six. The distance between two 
orthonormal bases vanishes when the bases are the same, and the distance reaches its maximal value 
of unity when the bases are unbiased. We perform a numerical search for the maximum average 
distance and find it to be strictly smaller than unity. This is strong evidence that no four mutually 
unbiased bases exist in dimension six. We also provide a two-parameter family of three bases which, 
together with the canonical basis, reach the numerically-found maximum of the average distance, 
and we conduct a detailed study of the structure of the extremal set of bases. 



I. INTRODUCTION 

Two orthonormal bases of a Hilbert space are said to 
be unbiased if the transition probability from any state of 
the first basis to any state of the second basis is indepen- 
dent of the two chosen states. In the finite dimensional 
case of C*, the normalization of the two basis states \ai) 
and I bj ) of two unbiased bases implies the defining prop- 
erty 



\{a,\b,)\^ = - foralH,j = l,2,...,d. 



(1) 



This maximum degree of incompatibility between two 
bases P, [1] states that the corresponding nondegener- 
ate observables are complementary. Indeed, the techni- 
cal formulation of Bohr's Principle of Complementarity 
Q that is given in Ref. Q relies on the unbiasedness of 
the pair of bases. Textbook discussions of this matter 
can be found in Refs. [5, 6], and Ref. [7] is a recent re- 
view on mutually unbiased bases (MUB), which are sets 
of bases that are pairwise unbiased. 

In addition to playing a central role in quantum kine- 
matics, we note that MUB are important for quantum 
state tomography [1, Q , for quantifying wave-particle du- 
ality in multi-path interferometers [lO|, and for various 
tasks in the area of quantum information, such as quan- 
tum key distribution [IL] or quantum teleportation and 
dense coding [T2| - [l^ . In the context of quantum state 
tomography, d -I- 1 von Neumann measurements provide 
d — 1 independent data each in the form of d probabili- 
ties with unit sum, so that in total one has the required 

— 1 real numbers that characterize the quantum state. 
A set of d -f 1 MUB is optimal, in a certain sense [^, for 
these measurements — if there is such a set. Such a set is 
termed maximal] there cannot be more than d -I- 1 MUB, 
since there are at most d + 1 (d — l)-dimensional orthog- 
onal subspaces in a (d^ — l)-dimensional real vector space 

Ivanovic [8| gave a first construction of maximal sets 
of MUB if the dimension d is a prime, and Wootters and 
Fields [9|] succeeded in constructing maximal sets when 
d is the power of a prime. These two cases have been 
rederived in various ways; see Refs. p"5l - [l7j . for example. 



For other finite values of d, maximal sets of MUB are 
unknown, but it is always possible to have at least three 
MUB (see 0] and references therein). 

The smallest non-prime-power dimension is d = 6. Lit- 
tle is known for sure about the six-dimensional case, for 
which Zauner has conjectured that no more than three 
MUB exist [31 . Numerical studies seem to support Za- 
uner's conjecture [H, [l^l- Computer-aided analytical 
methods, like Grobner bases or SemiDefinite Program- 
ming, have also been applied to this problem [2l[, but 
limitations in computational power have so far prevented 
any definitive answer. 

Recently, Bengtsson et al. 22] introduced a distance 
between two bases for a quantification of the notion of 
"unbiasedness." The distance vanishes when the two 
bases are identical and attains its maximal value of unity 
when they are unbiased. One can then consider the av- 
erage squared distance (ASD) between several bases and 
search for its maximal value. Importantly, this ASD is 
unity if the bases are pairwise unbiased, and only then. A 
numerical search for the maximum of the ASD between 
four bases in dimension six can be performed. Actually, a 
numerical study on essentially the same quantity was re- 
cently carried out by Butterley and Hall [23] ■ In terms of 
the ASD, they found the surprisingly large but strictly- 
less-than-one maximal value of 0.9983. This is strong 
evidence that no more than three MUB exist in dimen- 
sion six. However, the set of bases behind this maximum 
value is not reported in Ref. [23]. 

It is the objective of the present paper to close this gap. 
In Sec. [H] we review the notion of Bengtsson et al. for the 
distance between bases. We perform a numerical search 
for the maximum ASD between four bases in dimension 
six and report, in Sec. IIIIl our results which confirm the 
maximum found by Butterley and Hall. We then provide 
a two-parameter family of three bases which, together 
with the canonical basis, reaches the numerically-found 
maximum, for which we give a closed expression. We 
study this family in detail in Sec. llVl and conclude with 
a summary and outlook. Some matters of a technical 
nature are reported in the Appendix. 
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II. A DISTANCE BETWEEN BASES 

The main goal of this paper is twofold. First we nuraer- 
ically search for the maximum value of the ASD between 
four bases in dimension six and see that we cannot ob- 
tain four MUB. And second, we provide a two-parameter 
family of three bases which, together with the canonical 
basis, reaches the numerically-found maximum. 

Following Bengtsson et al. [Ill, we consider two or- 
thonormal bases of kets of C^, a — {|aj)} and h = {\bj)}, 
and quantify their squared distance by 

i,j=l ^ ' 

= d^I E l(a.|fc,)P(l-|(a.|?>,)p). (2) 

Clearly, this distance is symmetrical, Dah ~ ^ha and van- 
ishes when the bases are the same, that is: when the two 
sets of projectors {|ai)(ai|} and are identical; 

the maximal distance is unity. Dab < 1; and this maxi- 
mum is reached if the bases are unbiased, | (a^ p = 1 jd, 
and only then. 

In the original reasoning by Bengtsson et ai, Dab is 
actually the chordal Grassmanian distance of two planes 
in the {cP — l)-dimensional real vector space associated 
with traceless hermitian operators in the d-dimensional 
complex Hilbert space. One can also view Dab as the 
Hilbert-Schmidt distance between two rank-d statistical 
operators in (S) that are in one-to-one correspon- 
dence with the d-dimensional bases. Consult Appendix [51 
for this matter. 

For a set of k bases, we have the ASD between the 
k{k — l)/2 pairs of bases, given by 

Owing to the normalization, we have D^ < 1 with D^ = 1 
only if the k bases are pairwise unbiased. 

With this notion of distance at hand, we can numer- 
ically search for the maximum ASD between four bases 
in dimension six and see whether we obtain D^ = 1, or 
in other words, if we can find four MUB. This search is 
the subject matter of the next section. 



III. NUMERICAL STUDY 

Our numerical approach relies on the mapping between 
1-qudit operators and 2-qudit states established in Chap- 
ter 3 of 0. Plus we use the steepest-ascent algorithm to 
find the maximum ASD between four bases in dimension 
six. Details of the numerical method are presented in 
Appendix Our numerical results are reported below. 



TABLE 1: Rate of success and CPU time (in seconds) for the 
steepest-ascent search for the maximum ASD. The absolute 
maximum of = 1 is always reached for d -I- 1 bases in di- 
mensions d = 2, 3, 4, and 5. As the seven-dimensional case 
illustrates, the difficulty of finding the global maximum in- 
creases rapidly with the dimension because there are many 
local maxima at which the steepest-ascent search can get 
stuck. We have also looked for the largest ASD between four 
bases in dimensions two to seven. We could not find four 
MUB in dimensions two and six. The CPU time refers to a 
Intel®Core™2 Duo CPU E6550 processor at 2.33 GHz, sup- 
ported by 3.25 GB of RAM. 
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A similar numerical study was recently performed by 
Butterley and Hall [23| who minimized 1 — with the 
so-called Levenberg-Marquadt algorithm. Our approach 
confirms the extremal value they found, and we also ex- 
hibit the structure of the four bases that maximize D^ 
for d = 6. 

We have used our code not only in dimension d = 6 
but also for other d values as a mean of benchmarking. 
We have run our code 2,500 times for the dimensions two 
to five, 10,000 times for the dimension six and 300 times 
for the dimension seven, both for k = d + 1 bases and for 
four bases. Our results are summarized in Table [1] 

Only in two cases, the maximum ASD docs not reach 
the upper bound of = 1. They are the cases of four 
bases in dimension two and six. 

At most three MUB can be constructed in dimension 
two. Thus the maximum ASD between four bases has to 
be less than one. This example is interesting because it 
can be analytically solved. In M.^, the four bases corre- 
spond to the tetrahedron, where each edge represents a 
basis. 

Importantly, we have searched for the maximum ASD 
between four bases in dimension six. We have found the 
largest value to be f^max = 0.9983. In the search for 
the global maximum, we have also found a few other lo- 
cal maxima whose frequencies of occurrence are reported 
in Figure [1] These results are consistent with those re- 
ported by Butterley and Hall [23| . We find the same local 
and global maxima with very similar frequencies. This 
is as expected because we have generated the four ran- 
dom bases from which the search proceeds in the same 
way as Butterley and Hall, using the same dedicated 
Matlab command. The two numerical methods are dif- 
ferent, however. We use the steepest-ascent algorithm 
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FIG. 1: Histogram of the maximum values of the ASD found 
during a numerical search for 10,000 randomly chosen initial 
four bases. The search converges to one of the local max- 
ima in about 30% of all runs, and to the global maximum of 
^max = 0.9983 for the other 70% of initial bases. 



while they employ the Levenberg-Marquadt algorithm for 
a nonlinear least-squares optimization. 

Since we consider four bases, there are six pairs of 
bases and their respective distances are not without inter- 
est. Indeed, it turns out that one basis is unbiased with 
the three remaining bases. And these three remaining 
bases are themselves equidistant. The immediate impli- 
cation is that the privileged basis can be chosen to be the 
computational basis while the three remaining bases are 
Hadamard bases, that is: the unitary matrices composed 
of the columns that represent the basis kets with refer- 
ence to the computational basis are complex Hadamard 
matrices divided by -s/B. We recall here that a complex 
Hadamard matrix is a d-dimensional square matrix sat- 
isfying the two conditions of unimodularity and orthog- 
onality di] 



HW 



1 for i,j 
dt . 



1, 



(4) 



Therefore, the unitary matrix H/ y/d has matrix ele- 
ments that can be related to a pair of unbiased bases: 
{ai\bj) = Hij/Vd. 

In addition to maximizing D^, our code also returns 
the four bases for which the maximum is achieved. Af- 
ter a bit of polishing — the set of four bases is not unique, 
since global unitary transformations yield equivalent sets, 
and the order of kets in each basis is arbitrary — this al- 
lows us to seek for the structure hidden behind the max- 
imum ASD. In the next section we will present a two- 
parameter family of three bases. The two parameters 
are two phases while the three bases are three Hadamard 



bases. We study in detail the properties of this fam- 
ily and show that, for some definite values of the two 
parameters, these three bases together with the canoni- 
cal basis reach the numerically-found maximum ASD of 
0.9983. This definite structure of the optimal four bases 
is our main result, with a closed expression for -D^max as 
a most-welcome bonus; see Eq. below. 

Harking back to Table [TJ we note that the best set 
of seven bases in dimension six has an ASD of 0.9849, 
short of unity by a mere one-and-a-half percent. For 
all practical purposes — those of state tomography, say — 
these seven bases are marginally worse than the imagi- 
nary seven MUB that no one has managed to find. 



IV. THE TWO-PARAMETER FAMILY 

Following Karlsson '25'] , we express the two-parameter 
family in terms of 2 x 2 block matrices where each of 
the nine blocks is itself a complex Hadamard matrix. 
Such 2x2 block matrices are called -ff2-reducible. The 
two-parameter family contains three bases, the fourth 
basis being the canonical basis. We will see that these 
three Hadamard bases are equidistant, that their deter- 
minants are identical, and that they belong to the so- 
called Fourier transposed family Fq . Finally, we will 
show that together with the canonical basis they reach 
the numerically- found maximum of the ASD. 



A. Parametrization 

We begin by defining a few quantities. We will need 
the third root of unity w = exp(i27r/3) as well as the 
following 2x2 matrices: 



Z = 



1 


, x = 


X* 




1 1 




-1 




X 




1 -1 



and 



T 



1 Ljt^ 
1 



(5) 



where t = exp(i9t) and x — exp(i9x) are two phases. 
Let us notice that T and F2 are themselves Hadamard 
matrices. 

The Hadamard matrices for the three bases are given 

by 



Ml = 



a: 

iuj*tZX*^ 





M2 



1 

71 



X 



1 

71 



F2 F2 F2 
F2 F2 O-J * F2 
T oj*T ujT 



F2 F2 

T LoT u*T 
T uj*T ujT 



V6 



N2 



4 



and 



and 



X* 



1 

7E 






X3N3 . 






-itZX^ 



1 

71 



T 
F2 



to 



F2 F2 
ojT uj*T 

F2 U)F2 



(6) 



In the above parameterization, we have introduced the 
matrices Xi and Ni, i = 1, 2, 3, which we wih address as 
dephasing and central matrices, respectively. The deriva- 
tion of this parameterization is explained in Appendix IB] 
The next section is devoted to proving the three prop- 
erties earlier mentioned. Before we turn to the proofs, 
we wish to point out that an additional relation between 
the two phases x and t exists, 



cos(26>:,) 



sm 



(7) 



It reduces the two-parameter family to a single- 
parameter family. Of course, as a subfamily, it conserves 
all the fundamental properties of the two-parameter fam- 
ily. Furthermore it still reaches the maximum ASD. 



B. Properties 

1. Equidistance 

A significant property of the three proposed Hadamard 
matrices is their equidistance. The relevant terms that 
appear in the distance between the two bases Ma and 
Mb (i.e., |(ai|5j)|) are the elements of the product ma- 
trix M\M^^ (i.e., {ai\hj)) in absolute value. Therefore, if 
the three product matrices mIm^, M^M^ and M^M^ 
have equal coefficients in absolute value, then the three 
bases Mi, M2 and M^ are equidistant. This is exactly 
what happens here. Indeed, we have the following cyclic 
structure: 



mIm^ = - 



and 



ai 02 as 
03 fli a2 
a2 as ai 



MlM, 



6 



M^M^ 



Cl C2 C3 
C3 Cl C2 
C2 C3 Cl 



bi 62 &3 
63 61 62 

62 ^3 h 



(8) 



where, on the one hand, the 2x2 submatrices ai, 62 
and C3 have the same coefficients in absolute value and, 
on the other hand, 02, 03, 61, 63, ci and C2 have the 
same coefficients in absolute value. More precisely, these 
matrices have the following forms (where the symbol 
stands for swapping the two diagonal elements). First, 



ai 



a P 
-/?* a* 



fli 



C3 



iu)t* P iLU*ta 
iioVa* iujt* P 



(9) 



Second, 



02 



as 



Cl 



e a;*7* 



W7 — e 

-5* 7* 



it*5 iujtj 
iLj*t*-f* -iuj*t*e* 



bl = 0,2; 

h ^ 0,3, 

C2 = Cl . 



(10) 



The various coefficients in Eqs. ([9|) and (fTO)) can be ex- 
pressed in terms of the two angles 9^ and 9t, 



a = 4 cos 



(0,)(l-c^rsin(0,)) , 
-2iuj*t (cos{2e,,) - 2cos(6't - 27r/3) sin(6'^)^ 
-2a;* cos(6l^) (^uj* + 2t* sm{0^)j , 
-2it(^cos{2e^) - 2 cos(6lt) sin(6i^)^ , (11) 
-2iu*t* (^cos(26':r) - 2 cos(6't + 27r/3) sin(6l^)^ . 
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When Eq. ^ is fulfilled, we have e = ui*6* and a few 
simplifications arise. We obtain 

a3 = ujZa2Z , 63 = ujZbiZ , and C2 — —Zc\Z , (12) 

for example. 

2. Determinant 

A direct calculation shows that 

Det(A:i) = Det(iVi) = Det(X3) = Det(iV3) = wt^ . 

(13) 

Accordingly, the three Hadamard bases share the same 
determinant 

Det(Mi) = Det(M2) = Det(M3) = w*t'^ . (14) 

However, although the determinants are equal, there 
seems to be no simple relation between the three matri- 
ces Ml, M2, and M3. In particular, they do not have the 
same spectrum and are, therefore, not related by unitary 
operators. 



3. Fourier transposed family 

The Fourier transposed family, first studied by 
Haagerup ^26)], is parameterized by Karlsson in the 



5 



form |2£ 



F2 F2 F2 

Ti LoTi Lo*Ti 

T2 L0*T2 UJT2 



(15) 



where the 2x2 Hadamard matrices Ti and T2 are given 
by 



T,, = 



1 U 
1 -U 



(16) 



The equivalence relation in Eq. (jlSp means equality up to 
left and right dephasing and left and right permutations. 
In other words, the central matrix is the fundamental 
object that specifies the equivalence class. In the form of 
Eq. (|6]), it is clear that the three matrices A^i, N2, and 
N3 belong to the Fourier transposed family. As a result, 
the two-parameter family itself belongs to the Fourier 
transposed family. 

Let us note here that only the right equivalence is nat- 
ural for more than two bases as it states that bases are 
defined up to permutations and global phases of their 
basis states. In particular, the distance between bases 
is invariant under right equivalence but not under left 
equivalence. 




-3-2-10 1 2 3 

FIG. 2: Contour plot of the ASD for the two-parameter 
family. Along the dashed curves, relation ([7]) holds. The 
four single-parameter families — one for each dashed curve — 
are equivalent to the two-parameter family in the sense that 
the maximal and minimal value of the ASD can be found 
by searching along one of the dashed lines only. The ar- 
rows point to the location of one of the eight maxima at 
{Ox,et) = (0.9852, 1.0094), marked by a cross. 



C. Average distance 

Let us now compute the global maximum of the ASD 
between the three bases. Since the three bases are 
equidistant, we only have to compute the distance be- 
tween, say. Ml and M2. A direct calculation leads to the 
following expression 

Dl2i0^,et) = ^[5- P(sin(0,),cos(0t + iTT))] , (17) 

with the polynomial 

P{p, q) = 8/ -I- 8g^/ - IQq^p^ 

-f IQqp^ - 16gV + 8(7^/ 
- 7p^ - Uqp^ + Sq^p'^ 



+ 2p2 + Aqp . 



(18) 



We denote by (popt, ^opt) the {p, q) pair for which P{p, q) 
is minimal and, therefore, Di2{0x,9t) is maximal. It 
turns out that <7opt is related to popt by 



cos(0°P' + \tt) = go; 



1 - 2pgpt 

Popt 



(19) 



which is a particular evaluation of the function defined 
2 

opt 



in Eq. ^ , and p^p^ is the unique real solution of a cubic 
equation. 



112p^pt-192p^pt + lllPopt = 22, 



(20) 



that is 



sin(6'°' 



opt\2 



-Popt 



16r - 



28r 



0.6946 



(21) 



with r = (2lV3 - 36)^/''' = 0.7199. It follows that there 
are eight optimal pairs of phases (0°^*,^°^') for which 
the maximal distance D™2^ is reached. The above ex- 
pressions for 6'°P* and 6°^*" can be injected back into the 
formula of the distance to obtain first ^^2^ and then 



= — [71 - 12cos( 
70 ^ ^ 



1 

70 



71-3 



jopt)4] 

-12r-3 
"14^^ 



= 0.9983 , (22) 



which agrees with the numerically-found maximum ASD 
within the machine precision. 

Furthermore, the distance D12 vanishes for 



0. = 7r/2 , 9t 
and Ox = —it/2 , 6t 



(mod 27r/3) 

7r/3 (mod 27r/3) . (23) 



As can be verified from the parameterization ([6]) or from 
the matrix products ([8|), the bases are indeed identical 
up to global phases and permutations for these values of 
the two phases Ox and Ot- 

We can also consider the single-parameter family that 
we obtain when eliminating Ot by using Eq. ([7]). Since 
Eq. (fTO|) is equivalent to Eq. (O, this single-parameter 
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family reaches the maximum of the ASD — and also 
the minimum since Eq. ([T5|) is obeyed by {9x,0t) — 
(7r/2, 27r/3). This is illustrated in Fig. [21 a contour plot 
of D"^ for the two-parameter family of Hadamard bases, 
with the location of the {dx,Ot) values of the single- 
parameter family indicated. The location of one of the 
eight maxima is marked, and the locations of the other 
seven follow from the symmetry properties of the con- 
tours. 



V. SUMMARY AND OUTLOOK 

We performed a numerical search for the maximum 
ASD between four bases in dimension six. We found 
that it is strictly smaller than unity and so confirmed 
the recent study by Butterley and Hall [l^. We regard 
this result as strong evidence that no four MUB exist in 
dimension six. 

Next, we went beyond this numerical result by provid- 
ing the four bases behind the numerically-found maxi- 
mum. More specifically, we found a two-parameter family 
of three bases, which together with the canonical basis, 
reaches the maximum of the ASD. We characterized this 
two-parameter family in full. We proved its inclusion in 
the Fourier transposed family and shown that the three 
base are equidistant. Furthermore, we analytically com- 
puted the maximum ASD between these three Hadamard 
bases and the canonical basis to show that it reproduces 
the numerical result. 

Two directions might be relevant for an extension of 
the present study. First, it would be interesting to see if 
the optimality of our solution can be extended to a larger 
family of bases, for example, to the whole Fourier trans- 
posed family. Second and complementarily, there might 
exist an argument to restrict the search for the maximum 
ASD between the canonical basis and three Hadamard 
bases to the Fourier transposed family, instead of the en- 
tire Hadamard family which, so far, has not been fully 
parameterized. In this context, however, it should be 
noted that — as follows from the findings of Jaming et al. 
[27} — there are no four MUB if one restricts the search 
to members of the Fourier family. 

Finally, if there is no complete set of seven MUB in 
dimension six, the optimal measurement for state tomog- 
raphy, in terms of statistical errors, remains to be found. 
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Appendix A: Numerical method 

As discussed in Sec. 3.1 of Ref. 0, for any ket \ip) 
or bra in a d-dimensional Hilbert space H or 
respectively, there is a conjugate bra or ket 

9 (01 ^ 10*) e H (Al) 

such that 

= {^id^r = {^\^) . (A2) 

This mapping is not unique, but two different realizations 
differ at most by a unitary transformation. As a rule, {(f)* \ 
and (01 = \(j>)^ are different bras. 

Once a particular choice of mapping has been made, 
there is a one-to-one correspondence between one-qudit 
operators and two-qudit kets, 

\(p){<l>\eB{H)^\<l)*,^) en(E)n. (A3) 

In particular, for an orthonormal basis of kets inJi, a = 
{|ai), |a2) . . . , |ad)}, we have the conjugate basis a* = 
{ |a*), 102), . . . , jflrf) }, and jointly they are used in defining 
the two-qubit state 

1 

i=i 

which has the d-fold eigenvalue 1/d and the {(P — d)-iold 
eigenvalue zero. 

We normalize the Hilbert-Schmidt inner product of 
two-qudit operators in accordance with 

{A,B) = dTr{A''B} , (A5) 

so that {pajPa) = 1 and {pa,Pb) — ^/d for a pair of 
unbiased bases. For the two-qudit states associated with 
two single-qudit bases, we then have 

(P-P'') = ^ E \i^j\bk)\'^l~^Dl (A6) 

with the distance Dab of Eq. ([2]), where the identity 
I 1 2 

(0*0^1 6^. 5 J.) — \{aj\bk)\ is used. It follows that Dab can 
be expressed in terms of the Hilbcrt-Schmidt norm of 
Pa - Ph, 

Dab^^^^\\pa-Pb\\ (AT) 

with II A II = ^ {A, A). This tells us something important: 
If a 7^ fo, then pa 7^ Pb, so that the mapping a ^ pa is 
one-to-one. 

In passing, we note the following challenge. Clearly, 
not all two-qudit states with p = dp^ correspond to a 
single-qudit basis in the sense of Eq. (IA4p . But which 
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additional criteria identify the set of two-qudit states that 
do? 

We are interested in finding the maximum value of the 
ASD between k bases in dimension d. The numerical 
search begins with a randomly chosen initial set of bases, 
and then modifies the bases in each iteration round such 
that D'^ is systematically increased. 

An infinitesimal variation of a ket in basis a is given 

by 



S\aj) = iea\aj), 



where is an infinitesimal hermitian operator acting on 
the basis a. We have one such hermitian e operator for 
each basis. The resulting response of D"^ is 



<5Z?2^^tr{eaGa} 



(A9) 



a=l 



where tr{ } is a single-qudit trace and 



Ga 



(AlO) 

is the ath component of the gradient. If bases a and b 
are unbiased, there is no contribution to Ga from basis 
b and, therefore, there is no gradient for a set of MUB. 
But the converse is not true: We can have a vanishing 
gradient although the bases are not pairwise unbiased. 
When the gradient has nonzero components, we choose 
— nGa with a common k > that specifies the step 
size. This guarantees 5D'^ > if k is not too large, and 
maximization along the line specified by the direction of 
the gradient can be done by optimizing the value of k. 
The line optimization is a necessary ingredient if conju- 
gate gradients are used for accelerating the convergence; 
see Ref. [2S|, for instance. 



The finite unitary change of basis a, 
then accomplished by 



^ Vala," 



Va 

or K 



e 
1 



iea/2 



1 - ieal2 
or Va = {1 + tea) n [ 



1 + 



ri=0 



(All) 



or yet other ones, whichever of them is convenient to use. 
All three VaS equal 1 + iea to first order in Sa and differ 
in the higher-order terms. Note that a high-precision 
evaluation of the infinite product in the third version of 
Va requires very few terms. This makes the third version 
a viable alternative if the computation of the exponential 
in the first version or of the inverse operator in the second 
version is time consuming or imprecise. 

The iteration is terminated, when all components of 
the gradient vanish (in the numerical sense specified by 
the machine precision). We repeat this steepest-ascent 



search many times to ensure that we find the global max- 
imum. As Fig. [T] shows for {d,k) = (6,4), the iteration 
gets stuck in local maxima for about three attempts in 
ten and, see Table [U only four in ten trials are successful 
for (d,fc) = (6,7). 



Appendix B: Derivation of the two-parameter family 

The d X d matrix Uab composed of the transition am- 



(A8) 

plitudes {aj\bk) of two orthonormal bases is unitary. 



Uab = 



(«2| 

{ad\ 



|6i),|62),...,|&d) =Ul, 



UabUba = a • 



(Bl) 



The columns and the rows of Uab are representations of 
the kets \bk) and the bras (aj |, respectively. The unitary 
matrices associated with a set of bases have a compo- 
sition law for consecutive basis changes: Uab ~ UaJJcbi 
Uaa = 1. In particular, VdUab is a complex Hadamard 
matrix if the bases a and b are unbiased; see the para- 
graph containing Eq. (j?]). 

Now, from the numerical search we know that one of 
the bases that maximize the ASD between four bases in 
dimension six is unbiased with the other three bases. We 
identify this privileged basis as the canonical basis and 
refer to it as the zeroth basis, and characterize the set of 
four bases by the three 6x6 transition matrices 



Ml = Uoi , M2 = C/( 



02 : 



Ah = U 



03 : 



(B2) 



so that the columns of Mi are composed of the probability 
amplitudes of the kets of the ith basis with respect to the 
privileged basis. 

When multiplied by the matrices Mi, M2, and M3 
are 6x6 Hadamard matrices, for which we use Karlsson's 
parameterization ^25| . 

His parameterization applies to i?2 -reducible 
Hadamard matrices that can be written in the form 
H = XlPlNPhXj^, where the left and right X matrices 
only contain phases on the diagonal, the P matrices are 
permutation matrices, and the central matrix has the 
form 



N : 



F2 Ti T2 



(B3) 



where F2 is the unnormalized two-dimensional Fourier 
matrix of Eqs. ([5]) and the 2 x 2 matrices are those of 

Eq. (nr 



1 U 

1 -U 



with |ti| = 1 



(B4) 
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while 



(B5) 



with a unitary and hermitian 2x2 matrix A. It turns 
out that our Hadamard matrices are indeed iJ2 -reducible 
since they can be written as Mi — XLiPLiNiPn^XR. 
with the central matrices given by 



^1 = ^ 



1 

^7t 



1 

7! 



1 



F2 F2 F2 

F2 i^F2 6J*^2 

T Lo*T LoT 

F2 F2 F2 
T ujT uj*T 
T uj*T ujT 

F2 F2 F2 
T ujT lj*T 
F2 co"^ F2 UJF2 



(B6) 



see Eqs. 

As in Eqs. ([SJ, we choose to express the matrix T with 
factors oi oj — exp(z27r/3), 



T = 



1 ojt^ 
1 -ujt^ 



(B7) 



to exhibit the crucial dependence on the phase factor t. 
The left permutation matrices are all equal, P^^ = P^^ = 
Pl, = Pl- 

Third, we notice that only the left dephasing and per- 
mutation matrices are relevant for the distance. Indeed 
the right dephasing matrices only add global phases to 
the basis vectors while the right permutation only per- 
mute the basis vectors. In other words, two bases B and 
BPjiXji are equivalent in terms of distance. Therefore 
we can choose to conserve only the relevant structure for 
our bases, that is. Mi = XL-P^-Ni. 

The fourth step is to use the fact that only relative 
dephasing and permutations of the rows are relevant to 
the distance. Therefore we define new bases as 



Ml 
Mo 



Plxlx.P^N, 

N. 



2 ! 



M3 = Plxlx^P^N^ 



(B8) 



To simplify the notations, we again denote the two new 
diagonal matrices in P\^X\XiP]^ and Pj^X^X^P^ by Xi 
and X3, respectively. We further observe that 



Ai 
A2 
Ai 



and Xr, 



Si 
B2 
B3 



(B9) 



Next we add a suitable global phase to Xi and X3. 
We multiply Xi by exp(-iArg(Ai[l, l]Ai[2, 2]/2)) and 
X3 by exp(-iArg(Bi[l,l]Bi[2,2]/2)) such that Ai and 
Bi take the simple form 



exp(— 
ex.p{iq 



(BIO) 



for some phase (j). We end up with the remarkable form 



Xi = 



Ai 

^2 

Ai 



and X3 — 



AI 
uj*Al 
S3 



where [cf. Eqs. ([5])] 



X* 
X 



(Bll) 



(B12) 



So far, we have found that 



^3 

Bi — , 
B2 = uj*A 



^1, 
A* 



1 1 



(B13) 



and it only remains to find the structure behind the two 
2x2 dephasing matrices A2 and Bj,. 

To do so, we now consider the products Uij = mJMj. 
We obtain 



ai 


0,2 


as 




ai 




'F2AIF2 


0-3 


ai 


a2 


with 




= Fs 


F^A*T 


_a2 


as 


ai_ 




_a3_ 




T^AIT 



MIM2 



(B14) 

and is the standard (unnormalized) 3-dimensional 
Fourier matrix 



Fs 



1 1 1 

1 CJ CJ* 
1 LO* UJ 



(B15) 



Similarly we have 



and 





&2 63 




61 




F2B1F2 


bs 


61 62 


with 


62 


= Fs 


T^B2T 


b2 


63 bi_ 




bs 




pBsF2_ 



(B16) 



Cl 


C2 


C3 




Cl 




F2YIF2 


C3 


Cl 


C2 


with 


C2 


= Fs 


T^Y2F2 


_C2 


C3 


Cl_ 




_C3_ 




F2YST _ 



(B17) 
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where 





Yi 










A\ 








Y = x;x, = 





Y2 










UJAIA2 













^3. 













(B18) 

The seventh step is to look once more at the numerics. 
With respect to the product M-J^Mj, we see that 

a-2.=uj*ZaiZ . (B19) 

Thus we are lead to define the matrix equation 

Ex=a2-uJ*Za-iZ = (B20) 

This only represents a system of three equations since 
= £'i[2,2]. In the same manner, we have for 

tj 



(B21) 



and £^2(1, 1] = -^2(2, 2] so that, here too, only three equa- 
tions are relevant. Finally, for M^M^, we obtain 



E-i^c\ + ZC2Z = 



(B22) 



and, owing to (w* - l)i;3[l,2] = t{\ - iJ)E-i\l, 1], again 
only three equations are relevant. We should mention 
here that there are other interesting identities within the 
products M\M^ , such as &2 = [ai + a| + Z(ai - a|)Z]/2, 
but they are much more complicated to handle and will 
not be necessary to achieve our parameterization. 

The eighth step is to solve the above nine equations. 
We obtain 

Ey{\,\\ : tr{Ai} ^trjAg}, 

Ey{\, 2] : Ai - 2w*r2yl2 + ^*t*'^ Aj, = rl , 

i;i[2,l] : cj*r2Ai - 2w*r2A2 + A3 = r'l. (B23) 

From the numerics, we know that r = r' and thus Ai — 
A^, which we already found by looking at the dephasing 
matrix Xi . Note also that the expression of the complex 
number r is not required. Furthermore we find 



E2[lA 
E2%1] 



tr{Bi} ^utv{B2] , 

ijj*t*'^Bi + wBa - 2wr2A3 = si , 

Bi + r^Sa - 2ujt*^B3 = s'l . (B24) 



From the numerics, we know that s — s'{— r) and thus 
Bi = 0-1-62, which we already obtained by looking at the 
dephasing matrix X3 . The next three equations are much 
more interesting. Indeed we have 



E3[2,2] 



2tr{yi} - LJ*tT{Y2} - ujtT{Y3} = , 
2tr{yi} -ujt*hT{Y2} -uj*thT{Y3} =0, 



t*^Y2 - Y3 



From the numerics, we know that u = and the last 
equation reduces to 

Y3^t*^Y2. (B26) 



Since Y2 — UJA1A2 and Y3 = B^Ai, the above equation 
directly translates into 

B3 = u}*t^Al . (B27) 

This last relation can be inserted in i?3[l, 1] and E3[2, 2], 
which become identical and can be written as 



2tr{yi} - {uj* +Lut*^)tT{Y2} =0. (B28) 

This equation will soon become Eq. ([T]). 

A last hint from the numerics is needed. We actually 
notice that 



Y,Y2Y3 = -1. 



(B29) 



As Y3 = t*'^Y2, we arrive at t*'^YiY^ = ~t so that 
ujt*AlA2 = ±iU, where = 1, that is, U = 1 or 
U = Z since it has to be diagonal. With the help of the 
numerics, we conclude that 



and consequently 



A2 = iioHZAl'^ 



B3 = -itZAl 



(B30) 



(B31) 



The final parametrization of the dephasing matrices is 
therefore given by 



X, = 



X3 



Ai 
iuiHZAf 
^1 





AI 
oj*Al 
-itZAj 



(B32) 



which are ingredients in Eqs. ([5]). 

Let us finally come back to Eq. (jB28p . We can now 
substit ute Yi = Af and Y2 = {iuj*tZAf){ujAi) = itZAl 
in Eq. (|B28p and, upon defining x = exp{i0x) and t = 
exp(z0t), we arrive at 



(B25) 



cos(6't - 27r/3) 
which is Eq. ©. 



cos(26'j;) 



(B33) 
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